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ABSTRACT 

We study the effects of the sedimentation of the trace element ^^Ne in the coohng of white dwarfs. 
In contrast with previous studies — which adopted a simphfied treatment of the effects of ^^Ne 
sedimentation — this is done self-consistently for the first time, using an up-to-date stellar evolutionary 
code in which the diffusion equation is coupled with the full set of equations of stellar evolution. 
Due the large neutron excess of ^^Ne, this isotope rapidly sediments in the interior of the white 
dwarf. Although we explore a wide range of parameters, we find that using the most reasonable 
assumptions concerning the diffusion coefficient and the physical state of the white dwarf interior the 
delay introduced by the ensuing chemical differentation is minor for a typical 0.6 white dwarf. 
For more massive white dwarfs, say Mwd ^ 1.0 M©, the delay turns out to be considerably larger. 
These results are in qualitatively good accord with those obtained in previous studies, but we find 
that the magnitude of the delay introduced by ^^Ne sedimentation was underestimated by a factor of 
~ 2. We also perform a preliminary study of the impact of ^^Ne sedimentation on the white dwarf 
luminosity function. Finally, we hypothesize as well on the possibility of detecting the sedimentation 
of ^^Ne using pulsating white dwarfs in the appropriate effective temperature range with accurately 
determined rates of change of the observed periods. 

Subject headings: dense matter — diffusion — stars: abundances — stars: interiors — stars: evolution 
— stars: white dwarfs 
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1. INTRODUCTION 

White dwarf stars are the final evolutionary stage for 
the vast majority of stars and, hence, they play a key 
role in our quest for understanding the structure and 
history of our Galaxy. Standard stellar evolution theory 
predicts that most white dwarfs are the descendants of 
post-asymptotic gianch branch stars of low and interme- 
diate masses that reach the hot white dwarf stage with 
hydrogen-rich surface layers. The basic inner structure 
expected in a typical white dwarf consists of a degenerate 
core mostly composed of a mixture of carbon and oxygen 
(the ashes of core helium burning resulting from the evo- 
lution of the progenitor) plus some impurities resulting 
from the original metal content of the progenitor star. 
The most important of these impurities is ^^Ne, which 
results from helium burning on ^^N — built up during 
the CNO cycle of hydrogen burning — and reaches an 
abundance by mass of A"no ~ -^cno ~ 0.02 in Population 
I stars. 

An accurate determination of the rate at which white 
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dwarfs cool down constitutes a fundamental issue and 
provides a independent valuable cosmic clock to deter- 
mine ages of many Galactic populations, including the 
disk (Winget et al. 1987; Garcfa-Berro et al. 1988; Her- 
nanz et al. 1994) and globular and open clusters (Richer 
et al. 1997; Von Hippel & Gilmore 2000; Hansen et al. 
2002; Von Hippel et al. 2006). Thus, considerable effort 
has been devoted to observationally determine the lumi- 
nosity function of field white dwarfs — which also pro- 
vides a measure of the cooling rate — and to empirically 
determine the observed white dwarf cooling sequence in 
stellar clusters. Additionally, and from the theoretical 
point of view, the development of very detailed white 
dwarf evolutionary models, that incorporate a complete 
description of the main energy sources, has also been a 
priority, since it allows a meaningful comparison with the 
increasing wealth of observational data (Fontaine et al. 
2001; Salaris et al. 2000; Althaus & Benvenuto 1998). 

White dwarf evolution can be basically described as 
a simple cooling process (Mestel 1952) in which the de- 
crease in the thermal heat content of the ions consti- 
tutes the main source of luminosity. The release of both 
latent heat (Van Horn 1968; Lamb & Van Horn 1975) 
and of gravitational energy due to a change in chemi- 
cal composition during crystallization (Stevenson 1980; 
Garcfa-Berro et al. 1988; Segretain et al. 1994; Isern et 
al. 1991, 2000) also affects considerably the cooling of 
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white dwarfs. In particular, compositional separation at 
crystallization temporarily slows down the cooling. This, 
in turn, influences the position of the cut-off of the disk 
white dwarf luminosity function (Hernanz et al. 1994), 
which is essential in obtaining and independent determi- 
nation of the age of the Galactic disk. 

Another potential source of energy arc the minor 
species rich in neutrons, such as ^^Ne. Indeed, as first 
noted by Bravo et al. (1992) the two extra neutrons 
present in the ^^Ne nucleus (relative to A = 2Z) results 
in a net downward gravitational force and a slow, dif- 
fusive settling of ^^Nc in the liquid regions towards the 
center of the white dwarf. The role of ^^Ne sedimentation 
in the energetics of crystallizing white dwarf was first ad- 
dressed by Isern et al. (1991). The extent to which the 
cooling of a white dwarf could be modified by this pro- 
cess was later investigated by Bildstcn & Hall (2001) and 
quantitatively explored by Deloye & Bildsten (2002). De- 
pending on the underlying uncertainties in the physics of 
the interdiffusion coefficients, Deloye & Bildsten (2002) 
concluded that ^^Ne sedimentation could release suffi- 
cient energy to affect appreciably the cooling of massive 
white dwarfs, making them appear bright for very long 
periods of time, of the order of 10^ yr. Deloye & Bild- 
sten (2002) predicted that the possible impact of ^^Ne 
sedimentation on white dwarf cooling could be seen in 
metal-rich clusters. In fact, white dwarfs resulting from 
progenitors formed in metal-rich systems are expected 
to have larger abundances of ^^Ne in their cores, and 
the delay in the cooling history of white dwarfs resulting 
from ^^Ne diffusion would be largely amplified in these 
clusters. 

The recent detection of a markedly bright white dwarf 
population in the old open cluster NGC 6791 (Bedin et 
al. 2005) is certainly promising in this regard. Indeed, 
the significant supersolar metallicity characterizing NGC 
6791 (about [Fe/H] ^ +0.4) turns this cluster into a 
unique target to test the predictions of Deloye & Bild- 
sten (2002), who claim that the effects of sedimentation 
would be largest in this cluster. Although the anoma- 
lously bright white dwarf luminosity function observed 
in NGC 6791 could be reflecting the presence of a large 
population of massive white dwarfs with helium cores 
instead of carbon-oxygen cores, as suggested by Hansen 
(2005), ^^Ne diffusion constitutes a viable explanation for 
the morphology of the observed white dwarf luminosity 
function that cannot be discarded (Kalirai et al. 2007). 

In view of these considerations, a full and consistent 
treatment of white dwarf evolution is required to ascer- 
tain the validity of the predictions of Deloye & Bildsten 
(2002) about the role of ^^Ne diffusion in cooling white 
dwarfs, since their flndings were obtained using a sim- 
plified treatment of white dwarf evolution. In particu- 
lar, although Deloye & Bildsten (2002) carried out a full 
numerical calculation of the diffusion process, the white 
dwarf cooling sequences were computed using a series 
of static white dwarf mechanical configurations in which 
the thermal evolution of the white dwarf was decoupled 
and, moreover, the authors did not mention whether or 
not the mechanical structure of the white dwarf was re- 
computed self-consistently at each time step. Another 
crucial point is that the white dwarf models of Deloye 
& Bildsten (2002) were constructed using a fully degen- 
erate equation of state everywhere (with no inclusion of 



ion contributions in the structure calculation). This as- 
sumption is valid for the white dwarf interior, where the 
degeneracy is high, but not for the outer nondegenerate 
parts of the star. In addition, the treatment of Deloye 
& Bildsten (2002) required as an additional input a rela- 
tionship between the luminosity of the white dwarf and 
the temperature of the core for which they adopted 
the results of Althaus & Benvenuto (1998) which, ad- 
ditionally, was assumed to be isothermal. As a result, 
the heat sources and sinks could only be evaluated glob- 
ally and not locally as we do here. Consequently, these 
authors could not obtain absolute ages but, rather, rela- 
tive delays with respect to the absolute ages derived by 
Althaus & Benvenuto (1998). The aim of the present 
work is, precisely, to fill this gap. Specifically, we present 
a set of full evolutionary and self-consistent white dwarf 
cooling sequences, which include an accurate treatment 
of time-dependent ^^Ne diffusion and the physical state 
of the white dwarf interior. Special emphasis is given 
to the treatment of the luminosity equation, particularly 
regarding the contribution to the energy balance of the 
terms resiilting from changes in the profiles of chemical 
abundances. Calculations are followed down to very low 
surface luminosities. The plan of the paper is the follow- 
ing. Section 2 contains details about the input physcis 
of our white dwarf evolutionary code and the main as- 
pects of the diffusion treatment. In Sect. 3 we present 
our results. Finally, Sect. 4 is devoted to discussing and 
summarizing our findings. 

2. INPUT OF THE MODELS AND EVOLUTIONARY 

SEQUENCES 

2.1. General description of the code 

In this work we have computed the evolution of white 
dwarfs self-consistently, taking into account the changes 
in the abundance profiles induced by the gravitational 
settling of ^^Ne. That is, we have solved the full set of 
equations describing the structure and evolution of white 
dwarfs with the luminosity equation appropriately mod- 
ified to account for the energy released by the changes in 
the core abundances as a consequence of ^^Ne sedimen- 
tation. Consequently, this study constitutes a notable 
improvement with respect to the approach adopted by 
Deloye & Bildsten (2002) to compute the white dwarf 
cooling sequences, where all energy sources were treated 
as global quantities. 

The evolutionary code employed in this study is the 
one we used in our study of the evolution of massive 
white dwarf stars (Althaus et al. 2007), modified to in- 
corporate the effects of ^^Ne sedimentation in the in- 
terior of white dwarfs. In particular, the equation of 
state for the high-density regime is very detailed and it 
is described in Segretain et al. (1994). It accounts for 
all the important contributions for both the liquid and 
solid phases. The release of latent heat upon crystal- 
lization — which is assumed to occur when the ion cou- 
pling constant reaches r= 180 (Stringfellow et al. 1990; 
Chabrier 1993) — and neutrino emission rates have been 
included following Althaus et al. (2007). Radiative opac- 
ities are those from OPAL (Iglesias & Rogers 1996) com- 
plemented at low temperatures with the Alexander & 
Ferguson (1994) molecular opacities. Chemical redistri- 
bution due to phase separation upon crystallization has 
not been considered in this work. Particular attention 
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is given to the treatment of the very outer layers of our 
models at advanced stages of evolution. This is crucial 
in determining absolute ages. In particular, compression 
of these layers constitutes the main energy source of the 
white dwarf when it enters the Debye regime. In the 
present calculations, the internal solutions of the Henyey 
iteration have been considered up to a fitting mass frac- 
tion of « 10-"A'/^. 

The initial white dwarf configuration from which we 
started our calculations of the cooling sequences corre- 
spond to hot white dwarf structures that were obtained 
following the artificial evolutionary procedure described 
in Althaus et al (2005). Because the impact of ^^Ne dif- 
fusion on the cooling times is irrelevant for the hot white 
dwarf regime, the precise election of the initial condi- 
tions bears virtually no relevance for the purposes of the 
present work. We compute the evolution of white dwarf 
models with stellar masses of 0.6 and 1.06 M©. More 
massive white dwarfs are expected to have cores com- 
posed mainly of oxygen and neon, the result of carbon 
burning in semidegenerate conditions in the progenitor 
star (Ritossa et al. 1996; Garci'a-Berro et al. 1997). 
The core chemical profile of our models consists of a 
predominant chemical element (either carbon or oxygen) 
plus ^^Ne with an initially flat abundance throughout the 
core. The abundance of ^^Ne results from helium burning 
on via the reactions ^"^Nia, 7)i8F(/3+)i«0(a, -ff^Ne. 
In this work, we explore the evolution of white dwarfs 
resulting from Population I star progenitors, i.e., charac- 
terized by ~ ^CNO ~ 0.02. Although the theory of 
stellar evolution predicts the core of most white dwarfs to 
be composed essentially by a mixture of carbon and oxy- 
gen, we assume either carbon- or oxygen-dominated cores 
to be consistent with the one-component plasma under- 
lying assumption in the derivation of the diffusion coeffi- 
cient. The outer layer chemical stratification of all of our 
models consists of a pure hydrogen envelope of 10~^Af* 
(plus a small inner tail) overlying a helium-dominated 
shell and, below that, a buffer rich in carbon and oxy- 
gen. The initial chemical stratification corresponding to 
an oxygen-rich 1.06 Af© white dwarf is displayed in Fig. 
[TJ The shape of the outer layer chemical profile is given 
by element diffusion at low luminosities. Nevertheless, 
diffusion occurring in the outer layers was switched off in 
the present calculations. The evolutionary stages com- 
puted cover the luminosity range from log(L/LQ) « 
down to —5.3. 

For a proper treatment of white dwarf evolution in a 
self-consistent way with the changes in the core chemical 
composition induced by ^^Ne sedimentation, the stan- 
dard luminosity equation for evolving white dwarfs (ther- 
monuclear reactions are neglected) 
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has to be appropriately modified — see Isern et al. (1997) 
and also Kippenhahn et al. (1965). In Eq. ([T|), de- 
notes the energy per unit mass per second due to neutrino 
losses and the other quantities have their usual meaning 
(Kippenhahn & Weigert 1990). To this end, we write 
the energy equation in terms of changes in the internal 
energy per gram (u) and density (p) 
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Fig. 1. — The initial chemical abundance distribution for our 
1.06 Mq oxygen-core white dwarf. 
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and we assume the white dwarf interior to be made of two 
chemical elements with abundance by mass Xi and X2 
{Xi + X2 = 1), where Xi refers to the ^^Ne abundance. 
Now, 
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and with the help of the thermodynamic relation 

\9p)t.x.^ p' p' Urj,,x,' 
Eq. ([2]) can be rewritten as 



(4) 



dLr 
dM\- 



T f dP 



CvT 



du 



Xi. 



(5) 



Finally, by taking into account that p = p(P, T, Xi ) 
and the definitions 



6^- 



Eq. ^ can be written as 
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where A is given by 



4 



A 



( du 



PS / dp 



or alternatively 
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The second and third terms of Eq. ^ are the well- 
known contributions of the heat capacity and pressure 
changes to the local luminosity of the star. The fourth 
term represents the energy released by chemical abun- 
dance changes. Although this term is usually small for 
most stages of stellar evolution — as compared with the 
release of nuclear energy (Kippenhahn et al. 1965) — 
it will play, as it will be clear below, a major role in 
the cooling of white dwarfs with diffusively evolving core 
chemical compositions. For neutron-rich species, such 
as in the problem of ^^Ne diffusion in the core of white 
dwarfs, the derivative {du/dXi)p^T is dominated by the 
electronic contributions. Thus, A becomes negative and 
the last term in Eq. Q will be a source (sink) of energy 
in those regions where diffusion leads to an increase (de- 
crease) in the ^^Ne local abundance. In simpler terms, 
an increase (decrease) in the ^^Ne abundance — or in the 
molecular weight per electron /ig — will force the elec- 
tron gas to release (absorb) energy. Finally, in Eq. ^ 
we have included a term that represents the contribution 
of latent heat of crystallization (« k^T per ion, being fee 
the Boltzmann constant). 



2.2. Treatment of i 

The evolution of the core chemical abundance distribu- 
tion caused by ^^Ne diffusion has been fully accounted 
for. In presence of partial gradients, gravitational and 
electric forces, the diffusion velocities satisfy the set of 
equations (Burgers 1969) 
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Here, pi, pi, Ui, Zi and Wi mean, respectively, the partial 
pressure, mass density, number density, charge and dif- 
fusion velocity for species i and N means the number of 
ionic species plus electron. The unknown variables are 
Wi and the electric field E. The resistence coefhcients 
are denoted by Kij. For the white dwarf core we are 
mainly interested in the gravitational settling, thus we 
drop off the first term in Eq. The core of our white 
dwarf models consists of a predominant ion specie — ei- 
ther ^^O or -^^C plus ^^Ne (these two species are denoted 
by subindices 2 and 1, respectively). We assume the elec- 
tron to have zero mass. Since, pi = Ai rii nip, Eq. ([9]) 
can be written in the form 
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Al rii wi + A2 ^2 0^2 = (12) 

Here g is the local gravitational acceleration, mp is the 
proton mass, and Ai the atomic mass number. The last 
equation is the condition for no net mass flow relative 



to the center of mass. ^From the last set of equations 
we find that the diffusion velocities are (the effect of the 
induced electric field is considered) 
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where / = (1 -|- Aini/A2n2){niZi + n2^2)/"-i-^i'i2^2, 
and we have introduced the diffusion coefficient D as 

D = , . (15) 

[ni + na) if 12 

For 22Nc, Al = 22 and Zi = 10, while for either ^^O 
or ^^C, A2/Z2 = 2. Thus, we arrive at 
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The negative sign for the ^^Ne velocity reflects the fact 
that this ion settles towards the center of the white dwarf. 
This is expected because of the two excess neutrons of 
the ^^Ne nucleus (in comparison with A — 2Z). Note 
that if ^^Ne is assumed to be a trace element (tii « 0) 
then (711 -|- 7i2)/7ii7i2/ ~ Zi and we recover the diffusion 
velocity of Bildsten & HaU (2001) and Deloye & Bildsten 
(2002), namely wjqe = —2mpgD/kBT. 

For the diffusion coefficient D we follow the treatment 
adopted by Deloye & Bildsten (2002) for self-diffusion co- 
efficient in one component plasmas. In the liquid interior 
D becomes (in units of cm^/s) 
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where F is the Coulomb coupling constant. For the re- 
gions of the white dwarf that have crystallized, diffusion 
is expected to be no longer efficient due to the abrupt in- 
crease in viscosity expected in the solid phase. Thus, we 
set D = in the crystallized regions. To avoid numerical 
difficulties, the value of D is reduced from Dg to in the 
solid phase over a small range of F values. To assess the 
consequences of possible uncertainties in the values of D, 
we have performed evolutionary calculations also for the 
case D = 5Ds- Additionaly, we have explored the less 
likely situation in which white dwarfs experience a glassy 
transition at high F values instead of crystallizing. In a 
glassy state, diffusion is expected to continue. For this 
regime, viscosity increases considerably with the conse- 
quent reduction in the diffusion coefficient. Following 
also Deloye & Bildsten (2002), the diffusion coefficient 
in the glassy state is reduced as Dgiass = Dg/f* with 

/* = 1 + 5.4 X 10-7 F2-3666 + 2.16 X 10-27 

The evolution of the chemical abundance distribution 
in the core of our white dwarf models caused by diffu- 
sion is described in terms of a time-dependent, finite- 
difference scheme that solves elemental continuity equa- 
tions for the number densities rii. The equation govern- 
ing the evolution of the number density of species i in 
the space variable r is 
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Fig. 2. — Diffusively evofving internaf -^^Ne profffos for tfie 
0.6 Mq carbon-rich white dwarf sequence in terms of the mass frac- 
tion at selected evolutionary stages. The ages of the models are 
indicated in the legend. The numbers in brackets are the values of 
logTcff for the corresponding models. The upper (bottom) panel 
displays the situation in the case in which the diffusion coefficient 
is set to D = Ds {D = bDs). 
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supplemented with appropriate boundary conditions. 
Here, the diffusive flux caused by ion density gradients 
and the drifting flux due to gravitational settling are con- 
sidered. To solve this equation, we follow the method 
described in Iben & MacDonald (1985). In particular, 
we discretize Eq. as 
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where the superscript means that quantities are 
to be evaluated at the beginning of the time step 
At and the subscript k refers to the spatial index. 
A^fe = - rl_,^2)/^ and r^^./^ = (r^ + rl_^_,)/2. 

Note that oji and D are also averaged between adjacent 
grid points. Eq. ([20]) is linear in the unknowns rii. We 
follow the evolution of two chemical species in the core of 
the white dwarf: either ^^C or ^^O, and ^^Ne. In order to 
compute the dependence of the structure and evolution 
of our white dwarf models on the varying abundances 
self-consistently, the set of equations describing diffusion 
has been coupled to the evolutionary code. 

3. RESULTS 

We begin by examining Fig. [2] which displays the dif- 
fusively evolving ^^Ne profile in the core of the 0.6 
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Fig. 3. — Same as Fig. [2]but for the 1.06 Mq carbon-rich white 
dwarf models. 



carbon-rich white dwarf sequence as a function of the 
mass fraction. The situation at four selected stages of 
evolution is depicted. Clearly, diffusion appreciably mod- 
ifies the ^^Ne profile but only after long enough time has 
elapsed, causing a strong depletion of its abundance in 
the outer region of the core, where the diffusion timescale 
is comparable to the characteristic cooling time of the 
white dwarf. This is more apparent in the case of very 
efficient diffusion, as illustrated in the bottom panel of 
Fig. [21 which shows the evolving profiles when D = 5Ds 
is adopted. In this case, by the time the white dwarf has 
started the crystallization process (at log Toff = 3.753) 
diffusion has completely depleted ^^Ne in the outer- 
most region (20% by mass) of the white dwarf, and has 
markedly enhanced the abundance of ^^Ne in the central 
regions of the star. As crystallization proceeds, it leaves 
recognizable imprints in the inner ^^Ne profile. This is 
because diffusion is not operative in the crystallized re- 
gion, thus forcing ^^Ne to accumulate at the crystalliza- 
tion front. As the crystallization front moves outwards, 
an apparently irregular profile is left behind, as clearly 
seen in the figure. The apparent wave in the ^^Ne abun- 
dance in the crystallized core is a numerical artifact — 
see the discussion about this issue in Deloye & Bildsten 
(2002). This irregular profile, that arises from the use of 
finite timesteps in the evolutionary calculation and be- 
cause diffusion is switched-off at the liquid-solid bound- 
ary which occurs over a finite range of F values, bears 
no relevance for the calculation of the cooling ages of the 
white dwarfs. 

The rate at which ^^Ne diffuses downwards increases 
markedly with gravity. Thus, we expect a more rapid 
sedimentation and a faster depletion of ^^Ne in the outer 
layers of massive white dwarfs. This is borne out by 
Fig. [21 which illustrates the ^^Ne profile evolution for 
various 1.06 carbon-rich white dwarf models. The 
results for D — 5Ds, shown in the bottom panel of this 
figure, are worthy of comment. Indeed, while the white 
dwarf remains in a liquid state (at effective temperatures 
log Toff > 4.12), the ^^Nc abundance distribution will be 
strongly modified by diffusion during the evolution. Note 
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in the ^^Ne abundance — the term AXi in Eq. (jG} — for the same 
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of Fig. |3] The thin line denotes the heat capacity contribution for 
the model with logT^g = 4.154. Note that in the outer regions of 
the core, ^■^Ne diffusion constitutes an energy sink {AXi > 0). 



that by the time the white dwarf begins to crystaUize, 
^^Ne has diffused so deep into the core that no trace of 
this element is found in layers even as deep as 0.4 M© 
below the stellar surface. Note also the marked accumu- 
lation of ^^Ne towards the central regions. Shortly after, 
at log Toff « 4.02, about half the white dwarf mass is 
crystallized and the ^^Ne distribution will remain frozen 
with further cooling. Because massive white dwarfs crys- 
tallize earlier than less massive ones, ^^Ne diffusion will 
end at much higher effective temperatures when com- 
pared with less massive white dwarfs. This is a critical 
issue regarding the cooling ages of massive white dwarfs, 
as it will be discussed below. 

It is clear from these figures that diffusion substan- 
tially changes the core abundance distribution of ^■^Ne 
in the course of white dwarf evolution. These changes 
are in good agreement with those reported by Deloye & 
Bildsten (2002). We thus would expect that the con- 
tribution to the local energy budget of the white dwarf 
stemming from the last term in Eq. ([6]) could play a 
significant role. We show this in Fig. [J] for the same 
1.06 A/q carbon-rich stellar models shown in the upper 
panel of Fig. [3l The figure clearly illustrates the fact that 
^^Ne diffusion constitutes a local source or sink of en- 
ergy depending on whether the local abundance of ^^Ne 
increases or decreases. In fact, note that for those re- 
gions where diffusion has increased the local abundance 
of ^^Ne, AX22 < 0. Thus, this term constitutes an en- 
ergy source in such regions — in fact, the increase in the 
molecular weight forces the electron gas to release en- 
ergy. In the outermost part of the core however, where 
^^Ne depletion occurs, AXi > 0, thus implying a sink 
of energy in those regions. To compare with the main 
standard contribution to Lr, we include in Fig. [4] the 
run of CpT — second term in Eq. ^ — as a thin solid 
line. We do this for the model with logTcff = 4.154, 
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the fractional mass for the same 1.06 Mq carbon-rich white dwarf 
models shown in the upper panel of Fig. |3] The bottom panel 
shows the same quantity in the case in which ^^Nc diffusion is not 
considered. 



which corresponds to the model shown with dashed line. 
We have chosen this model because it corresponds to an 
evolutionary stage just before the onset of crystalliza- 
tion. For this model, the energy contribution from ^^Ne 
diffusion near the center of the white dwarf is similar to 
the heat capacity contribution — of about 0.003 erg/g/s. 
Finally, the dot-dashed and dotted curves correspond, re- 
spectively, to models for which 50 and 75% of their mass 
has already crystallized. The ^^Ne accumulation at the 
crystallization front leads to a large local energy contri- 
bution, restricted to a very narrow mass range. Because 
diffusion stops in the crystallized interior, the ^^Ne profile 
remains frozen there, and the term AXi vanishes. This 
fact limits the extent to which ^^Ne diffusion constitutes 
an energy source for the star. 

^From the above discussion, it is expected that the 
^^Ne profile evolution affects the fractional luminosity of 
the white dwarfs. This is exemplified in the upper panel 
of Fig. [5l which shows the fractional luminosity Lr in 
terms of the fractional mass, Mr for the same 1.06 M© 
carbon-rich white dwarf stellar models shown in Fig. |4l 
The bottom panel displays the behavior for the case in 
which ^^Ne diffusion is not considered. In both cases, 
the imprint of the latent heat release upon crystalliza- 
tion is apparent in the models with log Toff = 4.018 and 
3.942. There is clearly an important contribution of the 
^^Ne diffusion to the luminosity budget of massive white 
dwarfs. This contribution is notably enhanced when a 
larger diffusion coefhcient is considered, as it becomes 
clear from Fig. [SI which shows the same quantity in 
the case in which D = 5Ds is adopted. The models 
plotted in this figure are the same already shown in the 
bottom panel of Fig. \3\ Note that the fractional lumi- 
nosity profile is indeed strongly modified in this case. It 
is particularly relevant the model with logTcff = 4.115 
(dot-dashed line), which has already a small crystalline 
core. In fact, its internal luminosity at Air ~ 0.3M^ 
is about twice the luminosity emerging from the surface. 
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mass for the same 1.06 A/q carbon-rich white dwarf models shown 
in the bottom panel of Fig. \3\ in the case in which we adopt 
D = 5Ds. 



Note also the strong decrease in the fractional luminosity 
in the outer layers of this model, which results from the 
depletion of ^^Ne in these layers — see the bottom panel 
of Fig. [3] — with the consequent energy absorption. 

To get a deeper insight of the importance of ^^Ne dif- 
fusion into the global energetics during the whole white 
dwarf evolution, we show in Fig. [7] the resulting luminos- 
ity contribution obtained by integrating at each effective 
temperature the term —AXi in Eq. ^ throughout all 
the star. The result (expressed in solar units) is shown 
in terms of the white dwarf effective temperature for 
the 1.06 cooling sequences. The upper panel corre- 
sponds to white dwarf sequences with carbon-rich cores. 
To highlight the role of the core chemical composition, 
the bottom panel shows the sequences with oxygen-rich 
cores. Several aspects deserve comments. To begin with, 
note that, because the high photon luminosity of the star, 
the luminosity contribution from ^^Ne sedimentation is 
of a very minor importance during the hot white dwarf 
stages. It is only after the onset of core crystallization 
that this energy source contributes appreciably to the 
star luminosity. This is particularly true for the case of 
carbon-rich models, which cyrstallize at lower surface lu- 
minosities than their oxygen-rich counterparts. Because 
the rapid crystallization of the core, ^^Nc diffusion lumi- 
nosity declines steeply with further cooling. This decline 
occurs earlier in the case D = hDg because of the rapid 
depletion of ^^Ne in the outer layers in this case (see Fig. 
[3]). A similar behavior was reported by Deloye & Bild- 
sten (2002). As illustrated in the bottom panel of Fig. 
[3 the impact of ^^Ne settling in the energetics of white 
dwarfs with oxygen-rich cores is less relevant than in the 
case of carbon-rich ones. Finally, for both core compo- 
sitions, we compute the ^^Ne settling luminosity for the 
case in which the white dwarf is assumed to experience 
a transition to a glassy state at high F values instead 
of crystallizing. We show this for the cases D = -Dgiass 
and D = 5-Dgiass (dot-dashed and dot-dot-dashed lines, 
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The solid line displays the surface luminosity. The vertical line 
marks the approximate value of the effective temperature for the 
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respectively). Because diffusion is operative in a glassy 
state, ^^Ne settling luminosity obviously persists in more 
advanced stages of the evolution, providing the bulk of 
star luminosity at small effective temperatures. This is 
particularly true for the case D = 5Z)giass- However, at 
very low effective temperatures, the viscosity increase in 
the glassy state prevents diffusion from occurring, and 
the luminosity contribution eventually declines. 

For low-mass white dwarfs, the impact of ^^Ne set- 
tling is expectedly much less noticeable, albeit not negli- 
gible, particularly in the case of carbon-rich white dwarfs 
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Fig. 9. — Surface luminosity versus age for the 1.06 A'/0 white 
dwarf sequences with carbon- and oxygen-rich cores (upper and 
bottom panel, respectively). The solid line displays the cooling 
times for the case in which ■^^Ne diffusion is not considered. Dotted 
and dashed lines display the results for D = Ds and D = 5Ds, 
while the remaining curves correspond to the cases when models 
enter a glassy state. For the latter situation, we show results for 
D = -Dgiass and D = SDgiass- 

and efficient diffusion (see Fig. [8]). Note that because 
^^Ne sedimentation is a siower process in iow-mass wiiite 
dwarfs, it will take long for ^^Ne settling to contribute to 
the energy budget of the star. Since less massive white 
dwarfs crystallize at much lower effective temperatures 
than massive ones, this contribution will play some role 
at the very late stages of evolution. It is worth not- 
ing that the qualitative behavior of the evolution of the 
global energetics of our full sequences resembles that re- 
ported by Deloye & Bildsten (2002). In part, this is ex- 
pected since — as noted by Deloye & Bildsten (2002) 

— the diffusion timescale in the white dwarf core is 
much longer than the time on which energy is transferred 
through the core. We expect then a similar behavior be- 
tween both sets of calculations, as discussed below. 

The impact of ^^Ne sedimentation in the white dwarf 
cooling ages is seen in Figs. [9] and \W\ for the 1.06 
and 0.6 Mq white dwarf cooling sequences, respectively. 
Here, the white dwarf surface luminosity is shown as 
a function of the age. The solid line corresponds to 
the standard case when ^^Ne diffusion is not considered. 
The top and bottom panels of each of these figures are 
for carbon- and oxygen-rich cores, respectively. Clearly, 
^^Ne diffusion profoundly influences the cooling times, 
particularly those of massive white dwarfs. This influ- 
ence becomes more dramatic in the case of efficient dif- 
fusion — this is, the case in which D = 5Ds is adopted — 
which corresponds to the dashed curves. The star will 
spend a long time to get rid of the energy released by 
the diffusion-induced abundance changes, with the con- 
sequent marked lengthening of evolutionary times per- 
sisting until low luminosities. From the discussion in the 
preceding paragraph, it is evident that the signatures of 
^^Ne diffusion in the white dwarf cooling track start to 
manifest themselves earlier in more massive white dwarfs 

— at log(L/L0) « —2. The delay in the cooling times 
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Fig. 10. — Same as Fig. [9] but for the 0.6 M© white dwarf 
sequences. 

depends not only on the stellar mass, but also markedly 
on the chemical composition of the core. For instance, for 
the case D — 5Ds, the delays introduced in the 1.06 Mq 
carbon-rich sequence amount to 2.45 and 3.24 Gyr at 
log(L/L0) « —3 and —4, respectively, as compared with 
the 0.5 and 0.6 Gyr in the case of the oxygen-rich se- 
quence. In low-mass white dwarfs, appreciable delays in 
the cooling rates take place but only at low luminosities 
and for efficient diffusion. 

The cooling rate in the case that white dwarfs expe- 
rience a transition to a glassy state behaves differently 
since in this case ^^Ne sedimentation continues releasing 
energy even at very large F values. This distinct behav- 
ior can be appreciated in the two cooling curves of the 
1.06 Mq carbon-rich sequences shown as dashed and dot- 
dot-dashed curves which correspond, respectively, to a 
sequence that experiences crystallization (with D = bD^) 
and to a sequence that undergoes a glassy transition 
(with D = 5Dgiass)- In fact, the decline in the lumi- 
nosity contribution from ^^Ne sedimentation due to the 
presence of the crystalline core, see Fig. [71 translates into 
a change of slope in the cooling curve with D = 5Ds at 
log(L/L0) « —3. By contrast, the coohng rate for glassy 
white dwarfs remains much smaller. We thus expect sub- 
stantially larger delays in the cooling ages if a transition 
to a glassy state actually occurs in nature. In this case, 
massive white dwarfs would remain bright even at ex- 
ceedingly large ages. Note that during the crystallization 
stage, sequences that experience a glassy transition (in 
the case D = Ugiass) are younger than sequences without 
diffusion. This is because the glassy white dwarf does not 
release latent heat. Because crystallization is more rele- 
vant in less massive white dwarfs, this behavior is more 
evident in such white dwarfs. 

In closing, we note that the more sophisticated and 
self-consistent treatment of white dwarf evolution with 
diffusively evolving ^^Ne abundances done in this work 
yields results that are in qualitative agreement with those 
of Deloye & Bildsten (2002) which were obtained using 
a global and much more simplified treatment of white 
dwarf evolution. However, there are quantitative differ- 
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Fig. 11. — Difference in evolutionary times between sequences 
with various assumptions about diffusion and the sequence in which 
diffusion is neglected. Solid, dashed and dot-dashed lines corre- 
spond to sequences with D = 5Ds , D = Da , and glassy state with 
D = Dgiass ■ The thick (thin) set of lines are for the model sequence 
with carbon- (oxygen-) rich cores. 

ences. In the interest of comparison, we show in Fig. [TT] 
the age difference between our sequences with D = 5Ds 
and D = Dg and the sequence in which ^^Ne diffusion is 
not considered. Results for both carbon-rich and oxygen- 
rich core compositions are ihustrated. This figure can be 
compared with figure 10 of Deloye & Bildsten (2002). 
Although a direct comparison is not easy because Deloye 
& Bildsten (2002) consider cores composed by a mixture 
of carbon and oxygen, a close inspection of the figures re- 
veals that for both stellar masses our calculations predict 
age differences larger by a factor of ~ 2 when compared 
with the age differences obtained by Deloye & Bildsten 
(2002). In fact, we find that the impact of '^'^Ne diffu- 
sion on the white dwarf cooling ages is substantially more 
important than originally inferred by Deloye & Bildsten 
(2002). For the case in which a glassy transition is con- 
sidered, we find marked differences between both sets of 
calculations for the low-mass sequence. Indeed, in our 
calculations, the glassy sequence remains younger than 
the standard sequence without diffusion throughout the 
entire evolution (we stress that the glassy sequence does 
not release latent heat), while in the Deloye & Bild- 
sten (2002) treatment, once the crystallization process 
is finished, the released energy from ^^Ne diffusion in the 
glassy sequence yields an age increase. 

4. DISCUSSION AND CONCLUSIONS 

Motivated by the theoretical proposal that the cooling 
of white dwarfs could be altered by the sedimentation 
of ^^Ne in the liquid interior of these stars (Bravo 1992; 
Bildsten & HaU 2001; Deloye & Bildsten 2002) and the 
possibility that this could leave imprints in the white 
dwarf population of metal- rich systems, we have pre- 
sented in this paper detailed evolutionary calculations to 
address this issue. These are the first fully self-consistent 
evolutionary calculations of this effect. We have explored 
several possibilities. In particular we have thoroughly an- 



alyzed two classes of cooling sequences. The first class 
corresponds to an otherwise typical 0.6 white dwarf, 
whereas the second one corresponds to a massive 1.06 Mq 
white dwarf. Also the effects of the core composition 
have been addressed by following cooling sequences for 
carbon cores with a small admixture of ^^Ne and the 
corresponding ones in which the main constituent of the 
core was oxygen. The sensitivity of our results to the 
precise value of the rather uncertain diffusion coefficient 
has also been explored. Finally, we have also analyzed 
the effects of assuming that the white dwarf instead of 
crystallizing at low temperatures experiences a transi- 
tion to a glassy state. All in all, we have found that the 
sedimentation of ^^Ne has notable consequences in the 
evolutionary timescales of white dwarfs. The associated 
energy release notably delays the cooling of white dwarfs. 
The precise value of the delay depends on the adopted 
diffusion coefficient, on the mass of the white dwarf, on 
the core composition and on the possibility that white 
dwarfs experience a transition to a glassy state. The 
delay is larger for more massive white dwarfs and for 
carbon-rich cores. Specifically, for the case of an efficient 
diffusion the delay introduced in the 1.06 Mq carbon-rich 
sequence is 3.2 Gyr at log(L/L0) « —4 whereas it only 
amounts to 0.6 Gyr in the case of the oxygen-rich se- 
quence. In low-mass white dwarfs, appreciable delays in 
the cooling rates take place but only at much smaller lu- 
minosities. These results are in relative good agreement 
with the results of Deloye & Bildsten (2002), although we 
find that the delays obtained here are somewhat larger. 

To qualitatively assess the effect of ^^Ne sedimenta- 
tion on the white dwarf luminosity function, we have 
derived individual luminosity functions from our cool- 
ing curves. The number of white dwarfs of a given 
mass is proportional to the characteristic cooling time, 
n oc dt/d\og{L/LQ). We plot this quantity in Fig. [T2lfor 
the 1.06 and 0.6 Mq white dwarf cooling sequences with 
carbon-rich cores as a function of the surface luminos- 
ity. The characteristic cooling times in terms of the ef- 
fective temperature for the 1.06 white dwarf cooling 
sequences with carbon- and oxygen-rich cores are shown 
in Fig. [131 In both figures, the solid line corresponds 
to the case in which ^^Ne diffusion is not considered and 
the remaining curves to the cases of different diffusion as- 
sumptions, indicated in the legend. At large luminosities 
the signature of the accelerated evolution due to neutrino 
losses is clearly noticeable in the characteristic cooling 
times. At the extreme faint end of the luminosity range, 
the sharp decline in the luminosity function in the mas- 
sive sequences reflects the onset of Debye cooling, stage 
which is reached at different ages depending on the as- 
sumption about diffusion. 

Note that the effect of ^^Ne diffusion in the charac- 
teristic cooling times of low-mass white dwarfs with so- 
lar metallicity progenitors is barely noticeable. In sharp 
contrast, it is significative in the case of massive white 
dwarfs. In fact, depending on the internal composition 
and the uncertainties weighting upon the determination 
of the diffusion coefficient, an increase in the number 
of white dwarfs by a factor of « 5 could be expected 
in the effective temperature range where ^^Ne diffusion 
constitutes the main source of stellar luminosity. In the 
case that D = bDg and assuming carbon-rich composi- 
tion, this increase in the massive white dwarf population 
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would be expected at Tcs within 12,000 and 16,000 K. 
For oxygen- rich cores, only a modest increase is expected 
at somewhat larger effective temperatures. If a transition 
to a glassy state is assumed instead of crystallization, a 
marked increase in the luminosity function is expected in 
both carbon- and oxygen-rich massive white dwarfs. In 
the case of an oxygen-rich core, this increase takes place 
at the effective temperature range between « 12,000 and 
8,000 K. Interestingly enough, a possible large excess in 
the number of massive white dwarfs has been reported 
below Tcff = 12, 000 K (Kleinman et al. 2004; Liebert et 
al. 2005; Kepler et al. 2007). Although this increase is 
believed to be the result of a problem in the line fitting 
procedure, it might be possible that such an increase in 
the number of massive white dwarfs could in part be re- 
flecting a decrease in the cooling rate of massive white 
dwarfs, induced by the gravitational settling of ^^Ne in 
the core of these white dwarfs. 

The impact of ^^Ne sedimentation on the white dwarf 
structure can be assessed from pulsating white dwarfs. 
Indeed, as shown by Deloye & Bildsten (2002), the abun- 
dance gradient in ^^Ne caused by diffusion in the liq- 
uid regions leaves its signature in the Brunt- Vaisala fre- 
quency which, in turn, modifies the pulsation spectrum 
of ZZ Ceti stars at levels compared with the uncertainties 
of the measurements. More interestingly, pulsating white 
dwarfs with measured rates of period changes can be used 
to assess the impact of ^^Ne sedimentation on the white 
dwarf cooling directly. Indeed, the delay in the cooling 
times resulting from the ^^Ne diffusion is expected to al- 
ter the rate of change of the periods of ZZ Ceti stars. For 
a quantitative inference of the possible impact, we write 
the rate of change of the pulsation period of a pulsating 
white dwarf as 
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where T is the temperature of the isothermal core and a 
and h are constants of order unity which depend on the 
chemical composition, thicknesses of the hydrogen and 
helium envelopes, equation of state, and other ingredi- 
ents involved in the modeling of white dwarfs (Winget et 
al. 1983). For DA white dwarfs in the ZZ Ceti instability 
strip, the rate of change due to gravitational contraction, 
given by the second term of the right hand side of Eq. 
(|2T|) . is usually negligible and, thus, the secular rate of 
change of the period is directely related to the speed of 
cooling and it is a positive contribution. The rate of 
secular period change has been measured for some pul- 
sating white dwarfs. For instance, for G117-B15A, an 
intermediate mass ZZ Ceti star, Kepler et al. (2005) 
have derived the secular variation of the main observed 
period of 215.2 s, P = (3.57 ± 0.82) x IQ-^^ s s-\ with 
unprecedented accuracy. Other pulsating white dwarfs 
— like L 19-2 and R 548 — have also determinations of 
the secular rate of period change but are not as accu- 
rate as that of G117-B15A. We find that for the 1.06 Mq 
carbon-rich sequence with D = bDg the predicted period 
change decreases by a factor of about 1.7, as compared 
with the period change in the case that diffusion is not 
considered. This result is very preliminary, and an in 
depth study will be done in a forthcoming publication. 

In summary, although the effect of the sedimentation of 
^^Ne is sizeable its impact on the white dwarf luminosity 
function should be minor, except for a modest increase 
in the derived ages, for the most reasonable assumptions. 
However, the imprints of the sedimentation of ^■^Ne could 
be detectable using pulsating white dwarfs in the appro- 
priate effective temperature range with accurately deter- 
mined rates of change of the observed periods. 
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